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Abstract. We compute the small-tau expansion up to the third order for the form 
factor of two glued quantum star graphs with Neumann boundary conditions, by taking 
into account only the most backscattering orbits. We thus show that the glueing has 
no effect if the number of glueing edges is negligible compared to the number of edges 
of the graph, whereas it has an effect on the r 2 term when the numbers of glueing and 
non glueing edges are of the same order. 

PACS numbers: 05.45.Mt,03.65.Sq 
Submitted to: J. Phys. A: Math. Gen. 
1. Introduction 

The study of the Laplacian on a metric graph, a concept known as quantum graphs, now 
serves as a toy model for quantum chaos [H [2j [3]. Indeed, there exists an exact trace 
formula relating eigenvalues and periodic orbits. Moreover, depending on the graphs, 
exact computations of these orbits may be possible, whereas they are out of reach in 
most dynamical systems. It has thus been shown [U HI E] that spectral statistics of 
simple generic graphs follow random matrix statistics when the size of the graph tends 
to infinity, as expected for chaotic quantum systems. \^-star graphs with Neumann 
boundary conditions (graphs formed by a central vertex connected to V other vertices 
by edges of different lengths) play a special role because of the high degeneracy of 
their periodic orbits. As could be expected, this degeneracy breaks the random matrix 
statistics : this has been shown by the computation of the two-point correlation function 
([HI [?])• We will investigate here what happens when two star graphs are glued together. 
We will first show that glueing two star graphs with o(V) incommensurate bonds has no 
effect on the form factor when V tends to infinity, and we will next compute the small 
t expansion of two V^-star graphs glued by 0(V) bonds. 

Spectral statistics for quantum graphs and the trace formula relating them to 
periodic orbits will be presented in the first part, as well as a presentation of the model 
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we will be dealing with. The second part will recall the small r expansion of the form 
factor for star graphs as obtained in [7j. In a third part, we will explain our computations 
in the glued case. 

2. Quantum graphs : the trace formula, and the quasar model 

We will start by some vocabulary and notations. Let G = (E, V) be a graph with 
a metric structure : to each edge G E C V x V is assigned a length i^-, 

such that lij = Iji, although the graph is supposed to be non oriented, that is 
G E => G E, and kj = lj iy we will consider the edges to be oriented : 

is different form it really describes the edge going from i to j). On each 

edge one can thus define a coordinate x such that x = corresponds to the 

vertex i, and x = kj corresponds to the vertex j. A periodic orbit of period n is a set 
of n edges (pi, ■ ■ ■ ,p n ) such that pi ends where pi + \ starts (as well as p n and p±). A 
periodic orbit is called primitive if it is not the repetition of a shorter periodic orbit. 
A primitive orbit repeated r times is a non-primitive orbit with repetition number r. 
{J} will denote the equivalence class of all orbits of length I. On each edge one 
is looking for the spectrum of the Laplacian. In other words, one wants to find A and 
ipij such that —-^r- = \ 2 ipij(x). As one looks for eigenfunctions defined on the whole 
graph, one imposes continuity relations at each vertex, ipijifS) = ipik(0). Moreover, the 
function should have a unique value on a given point, regardless of the sense of the edge 
it belongs to : hence one wants ipij(x) = 4>ji{hj — x). This actually corresponds to the 
Neumann condition on each vertex . = 0. It is then a simple exercise to check 
that the eigenvalues A are the solutions of det(J — e~ tXL S) = 0, where S and L are 
\E\ x \E\ matrices : L is diagonal with the length of each edge as diagonal element, and 
S is defined by SVij),(j,fc) = —&i,k + rjr, where Vj is the number of j's neighbours. 

The trace formula as derived in |T] states that if d(X) = ^ n 5(X — A n ) is the spectral 
density, then d(X) = ^ + ^ ^nS P gp n 7"A> cos (A^). Here L is the total length of all 
edges, P n is the set of all periodic orbits of period n up to cyclic reordering (that is 
VqiV\iVi an d Pi,P2,Po are the same orbits), l p is the length of the orbit, r p its repetition 
number, and A p = f]" = i S PuPi+1 . 

A V-star graph is a graph with vertices {0, ..V} and edges E = {(0, i), (i,0), 1 < 
% < V} : the V vertices are all connected to the center 0. The S- matrix elements for 
such a graph are £(0,1), (i,o) — 1 (this corresponds to trivial scattering), S^o^cm) — — 1 + p= 
(backscattering) and S(i t o),(oj) = y (normal scattering). The lengths will be taken so 
that they are incommensurate, and that their distribution is peaked around 1 : for 
instance, they can be chosen randomly, uniformly in [1 — ^7, 1 + ^77], each length being 
independent from the others. 

We will consider two such Vi-star graphs and V^-star graphs, and connect them 
with M edges linking the two centers. We will call the resulting graph a (Vi,V2, M) 
quasar graph, and denote V — V\ + V 2 + M. The edges of the two star graphs will be 
denoted by roman letters (a, b, and so on) and the glueing edges by greek letters (e, (, 
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Figure 1. A quasar graph with V\ = V2 = 6 and M = 3 



and so on). 

The M glueing bonds are also supposed to have incommensurate lengths, chosen 
randomly in [1 — ^7,1 + ^y]- Hence a periodic orbit of period 2n has a length in 
[2n — y , 2n + y] ; such intervals for different n less than V do not overlap. Of course, in 
a periodic orbit of period 2n, each non glueing edge is visited an even number of times 
(and there is are trivial backscatterings) , whereas a glueing edge can be visited an odd 
number of times (but of course the total number of visits oof glueing edges has to be 
even) . The contribution A p of such an orbit will depend on the number of backscatterings 
and normal scatterings, and will also depend on the center on which the scattering takes 
place : it will be a product of factors r\ = (— 1 + Vi + M ) (f° r backscatterings on graph i) 
and ti = yr+M (f° r norma l scatterings on graph i). 

The interesting limit will be the limit V, n tends to infinity with r = y fixed. In 
this limit, orbits with the biggest contribution A p will be orbits with the largest number 
of backscatterings. We will suppose that V\ = v{V and V 2 = v 2 V , and consider two 
cases : M = u 3 V , or M = o(V). 

3. Form factor for the V-star graph 

The 2-point correlation function is defined as R 2 {x) = (^) 2 (d(\)d(\ - ^)). 
The brackets denote a mean value with respect to the A's, that is (/) = 
liniA^oo 7^- J^ A f(X)d\. Using the trace formula and performing the integral, one gets 



where the sum is over the pairs of periodic orbits (p,p'), up to cyclic permutations, of 
lengths l p and l p > and repetition numbers r p and ry. 

The form factor is defined by K st (r) = j^(R 2 (x) — 1) exp(2?7rxr)dr. When V is 




Small-tail expansion for the form factor of glued quantum star graphs. 



4 



large, K*(t) for r G [y — ^U, y + ^77] is well approximated by 




where the first sum is over isometry classes of periodic orbits of period 2n and the second 
is over a given isometry class. 

A combinatorial analysis of periodic orbits leads Q7J) to a formula valid near r = 0, 
yielding a small r expansion K st {r) = 1 — 4r + 8r 2 — |r 3 + o(r 4 ). 
One should note that the first three terms are given by the orbits consisting of only one 
edge; to compute the r 3 and r 4 terms, one also needs to take into account the isometry 
classes of orbits consisting of two or three different edges, but one can then only take into 
account in these isometry classes the orbits with the maximum number of backscattering 
(that is, orbits of the form aa...abb...b, or aa...abb...bcc...c). This approach starts only 
to fail for the r 5 term, whereas the diagonal approximation fails at r 4 ([6]). 

4. The form factor for the quasar graphs if 1/3 — 

We thus want to compute for r G [y — ^7, y + 7^7] fixed, 



K{t) = 




Some terms in the sum over equivalence classes concern periodic orbits visiting only 
one star graph (and no glueing edge). This yields two subsums very similar to the form 
factor of one star graph : the difference is that in the subsum over the i-th stargraph, the 
number of available vertices is Vi, whereas the scattering and backscattering factors are 

U = v-+m anc ^ r * = ( — 1 + v-+m )' anc ^ no ^ V anc ^ ( — 1 + f)" Nevertheless, if M — o(V), 
these two sums become in the limit (n, V) tends to infinity viK st {^) +z/ 2 X^(^-), where 
K st is the form factor for a star graph presented in the preceding section. 

We will compute the corrections of the most backscattering orbits to these unglued 
terms, and check that they disappear in the limit. Since these orbits should give the 
correct r — >• asymptotics, this means that the form factor of the quasar graph is in 
this case V\K st {^) + V2K st (^-), and that the glueing has no effect. In particular, if 
v\ = v-i = \, one gets K st (2r). 

Let us consider orbits of length around 2n with maximal backscattering, that is 
orbits consisting only of one glueing edge. Each equivalence class consists of one orbit 
a. ..a with repetition number n, and A p = (r 1 r 2 ) n . There are M such equivalence classes. 
Hence their contribution to K is -^M(2n) 2 ^ ri ^ which, when V, n tend to infinity such 
that y is finite, is of the order y-, hence negligible if M — o(V). This is not surprising 
: there are not enough such terms for them to contribute. This will be the case for all 
orbits containing glueing edge : to make up for the y terms corresponding to normal 
scatterings, one needs a V combinatorial factor, that one does not get when using glueing 
edges. 
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5. The z/ 3 7^ case. 

The situation is very different here, since even the two subsums do not give a proper 
star graph limit. We will thus here compute the exact contributions of the most 
backscattering orbits of length around 2n. 



5.1. Orbits with no scattering 

Let us first examine the case of the no scattering orbits. They consist of one edge, and 
they can be of three kinds, depending if the edge is a glueing edge or not. In all cases, 
the equivalence class is a singleton. 

If the edge is on star graph i, one gets A p = r™ and a repetition number n. There 
are Vi such equivalence classes, yielding in the limit Vi exp(— 4 ^ ) = z/jexp(— 4rj) if 
one writes = T . . 

If the edge is a glueing edge, A p = (rir 2 ) n (remember that there is no trivial 
backscattering for glueing edges). There are M such classes, giving a total contribution 
of z/ 3 exp(-4r(^ + ^)) = z/ 3 exp(-4(r 1 + r 2 )). 

5.2. Orbits with two scatterings 

These are orbits with two different edges. There are three different cases. 

If these two edges belong to the same star graph i, the weight A p is t 2 r™ -2 . There 
are Vi ^~^ ways to choose the edges, and n — 1 ways to choose how many times each 
orbit is visited. 

This yields when V tends to infinity a contribution to the form factor 

If the two edges are glueing edges, there are MiM — 1)/2 ways to choose them. But 
the weights depend on the parity of the number of visits of each edge. If it is odd, there 
is one scattering on star graph 1 and one scattering on star graph 2, hence a weight 
A p = ^^(^i^)" -1 - If h is even, the scatterings can happen either on star graph 1 or 
on star graph 2 : there are two orbits giving the main contribution in the equivalence 
class, with respective weights A p = t\r™~ 2 rV; and A p = t\r2~ 2 r™. Even or odd, there are 
n or n — 1 ways to choose how many times each orbit is visited. 

All in all, this case yields when V tends to infinity a contribution to the 
form factor of 8rV 3 2 exp(-4r 2 ) exp(-4r 1 )(^^ T ^ F + {j^i + j^^i?) = 

Last, if one of the two edges is a glueing edge and the other one is a i-star graph 

edge, there are MVi ways to choose the edges; there are n — 1 ways to choose the number 

of visits of each edge, but the weight depends on it : for example, if % — 1, and if the 

— 

glueing edge is visited n e times (n e has to be even), the weight is A p = tlr±~ r 2 2 . The 
sum over equivalence classes involves a sum over n e which is a geometric sum; in the 
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Figure 2. K*(2t) (circles), and the form factors for a (50,50,1) quasar graph 
(squares) and for a (50,50,50) quasar graph (stars). The solid line corresponds to 
the r 3 expansion in the z/3 = 0, vi — V2 case, the dashed line corresponds to the r 3 
expansion in the v\ = v-x = v% case. 



limit V — > 00 one gets a term 
_ 2 4z/ 3 z/i(z/ 2 + u 3 ) 



T 



exp(-4n)(l -exp(-4r 2 )) 



I6Z/3Z/1 



All in all, the first contribution of the two scatterings orbits to the small r expansion 
of the form factor is a r 3 term equal to 



+ 



1 



+ 



31/g 



Oi + z^) 2 (z^ + ^i) 2 {yi + ^) 2 (^2 + V3) 2 , 

5.3. r 3 expansion for the z/ 3 7^ case 

Putting everything together, the r 3 expansion of form factor for the quasar graph is 
l + 3z/ 3 § ^ + 14^^ + Uu 3 u 2 + 17 u 3 2 + u 2 2 



ir + i 



-T 



In the particular case v\ = is 2 , (and hence v 3 = 1 — 2v>i), one gets 

14i/f - 40z/i + 17 



-r 



1 



16(1-31/!) 2 
(l-^i) 2 3 



r 3 + o(t 3 ) 



The t 2 coefficient is then comprised between 32 (corresponding to z/3 = 0) and 36 
(corresponding to V\ — v% — z/3). The difference of the r 3 term with the unglued z/3 = 
case is also maximal for v\ — v% + 1/3. Yet, even in this case, the two expansions in the 
glued case and the unglued case are very similar for r £ [0, 1], so that it is difficult to 
give numerical evidence. As shown on figure [2j the two form factors are really close to 
each other until r ~ 0.15, and one would then need higher order terms. 
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6. Conclusion 

The first effect of the glueing is in the r 2 term : the non-universal statistics of the 
quantum graph seems to be quite robust. This is not so surprising, since this non 
universal statistics is strongly linked to the fact that the random walk on the graph 
converges very slowly to the uniform distribution, and this fact cannot be changed 
by the glueing. It should be expected that any "star-graph" component of a quantum 
graph, as long as it is macroscopically big (containing a non- zero fraction of the vertices), 
should yield a form factor such that K(0) = 1. 

This work benefitted from the support of the french ANR project "Arithmatrix" . 
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